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[57] ABSTRACT 

A digital image compression preprocessor for use in a 
discrete cosine transform-based digital image compression 
device is provided. The preprocessor includes a gathering 
mechanism for determining discrete cosine transform sta- 
tistics from input digital image data. A computing mecha- 
nism is operatively coupled to the gathering mechanism to 
calculate a image distortion array and a rate of image 
compression array based upon the discrete cosine transform 
statistics for each possible quantization value. A dynamic 
programming mechanism is operatively coupled to the com- 
puting mechanism to optimize the rate of image compres- 
sion array against the image distortion array such that a 
rate-distortion-optimal quantization table is derived. In 
addition, a discrete cosine transform-based digital image 
compression device and a discrete cosine transform-based 
digital image compression and decompression system are 
provided. Also, a method for generating a rate-distortion- 
optimal quantization table, using discrete cosine transform- 
based digital image compression, and operating a discrete 
cosine transform-based digital image compression and 
decompression system are provided. 

44 Claims, 5 Drawing Sheets 
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IMAGE COMPRESSION SYSTEM AND 
METHOD HAVING OPTIMIZED 
QUANTIZATION TABLES 


This invention was made with United States Government 5 
support awarded by the NASA. Grant No. NAGW-3914 and 
NSF, Grant No. IRI-9224741. The United States Govern- 
ment has certain rights in this invention. 

This invention was made with United States Government 
support awarded by the National Science Foundation (NSF). to 
Grant No. IRI-9224741 and the National Aeronautical Space 
Administration (NASA). Grant No. NAGW-3914. The 
United States Government has certain rights in this inven- 
tion. 


FIELD OF THE INVENTION 


15 


The present invention relates to digital image compres- 
sion and more particularly to compression techniques which 
utilize discrete cosine transforms. 

BACKGROUND OF THE INVENTION 


20 


The Discrete Cosine Transform (DCT) lies at the heart of 
many commonly used lossy image and video compression 
schemes such as the Joint Photographic Experts Group 
(JPEG). Motion Picture Experts Group (MPEG), and H.261 25 
compression standards. It will be appreciated by those 
skilled in the art that MPEG currently consists of several 
different algorithms including: the original MPEG algorithm 
which was designed for video (motion) images, the MPEG-2 
Video algorithm which was designed to specify a coded 30 
bit-stream for high-quality digital video, the MPEG-2 Audio 
algorithm which was designed for multi-channel audio 
coding, and the MPEG-2 Systems algorithm which specifies 
how to combine multiple audio, video, and private-data 
streams into a single multiplexed stream. ° 

The extent of compression achieved in any of these 
compression algorithms depends upon the coarseness of 
quantization of the transform coefficients. The coarser the 
quantization, the lesser the entropy of the quantized coeffi- 
dents. But. coarse quantization also leads to poor quality of 
the decompressed image. Thus, the particular quantization 
table used directly determines the rate-distortion tradeoff 
(i.e. . the compression-quality tradeoff). 

Several approaches have been tried in order to design 4J 
quantization tables for particular distortion or rate specifi- 
cations. The most common of these is to use a default table 
and scale it up or down by a scalar multiplier to vary quality 
and compression. Other approaches include psycho-visual 
model based quantization, rate-distortion model based 50 
quantization, and stochastic optimization techniques. 

A need still exists for a compression algorithm for opti- 
mum quantization table design that does not rely on visual 
or rate-distortion models and is not very expensive in terms 
of computation time. The algorithm should admit a wide 55 
range of quality measures such as Peak Signal to Noise Ratio 
(PSNR) and weighted PSNR. Also, the algorithm should 
produce quantization tables that optimize the tradeoff 
between quality and compressed size. In addition, the algo- 
rithm should simultaneously optimize quantization tables 50 
over a wide range, rates as well as distortions. 

The present invention provides a solution to this and other 
problems, and offers other advantages over the prior art. 


SUMMARY OF THE INVENTION 

A digital image compression preprocessor for use in a 
discrete cosine transform-based digital image compression 


device is provided. The preprocessor includes a gathering 
mechanism for determining discrete cosine transform sta- 
tistics from input digital image data. A computing mecha- 
nism is operatively coupled to the gathering mechanism to 
calculate an image distortion array and a rate of image 
compression array based upon the discrete cosine transform 
statistics for each possible quantization value. A dynamic 
programming mechanism is operatively coupled to the com- 
puting mechanism to optimize the rate of image compres- 
sion array against the image distortion array such that a 
rate-distortion-optimal quantization table is derived. In 
addition, a discrete cosine transform-based digital image 
compression device and a discrete cosine transform-based 
digital image compression and decompression system are 
provided. Also, a method for generating a rate-distortion- 
optimal quantization table, using discrete cosine transform- 
based digital image compression, and operating a discrete 
cosine transform-based digital image compression and 
decompression system are provided. 

These and various other features as well as advantages 
which characterize the present invention will be apparent 
upon reading of the following detailed description and 
review of the associated drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram showing a preferred embodi- 
ment discrete cosine transform-based digital image com- 
pression and decompression system in accordance with the 
present invention. 

FIG. 2 is a flowchart of the preferred embodiment 
RD-OPT algorithm used in the system shown in FIG. 1 in 
accordance with the present invention. 

FIG. 3 is a more detailed flowchart of the preferred 
embodiment DCT statistics gathering step in accordance 
with the present invention for the RD-OPT algorithm shown 
in FIG. 2. 

FIG. 4 is a more detailed flowchart of the preferred 
embodiment R„(q) calculating step in accordance with the 
present invention for the RD-OPT algorithm shown in FIG. 
2 . 

FIG. 5 is a more detailed flowchart of the preferred 
embodiment D„(q) calculating step in accordance with the 
present invention for the RD-OPT algorithm shown in FIG. 
2 . 

FIG. 6 is a more detailed flowchart of the preferred 
embodiment finding RD-optimal Q step in accordance with 
the present invention for the RD-OPT algorithm shown in 
FIG. 2. 


DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

The human visual system is not very sensitive to sudden 
changes in intensity across an image. Lossy image compres- 
sion techniques strive to discard that part of the image 
structure which is less perceptible to the eye. The two- 
dimensional Discrete Cosine Transform (DCT) offers an 
efficient way to break up the underlying structure of an 
image into different spatial-frequency components. The 
high-frequency components are less perceptible to the eye 
compared to the low-frequency components. Thus, the DCT 
orders the information-content of an image into parts with 
different degrees of visual importance. These parts can then 
be selectively discarded or stored, with varying degrees of 
accuracy, for lossy compression of the image. The 
compression-ratio increases as more and more information 
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is thrown away. DCT-based compression techniques typi- 
cally allow the user to specify a table (called the quantization 
table) that stipulates the level of accuracy with which each 
spatial-frequency component is to be stored. 

5 

The basic steps used in DCT-based image compression 
and decompression will be described below. 

Let I be a WxH image with pixel values in the range [0 
. . . MJ. The DCT-based compression process consists of the 
following steps. to 

1. The image is divided into 8x8 blocks. To each image 
block f, the DCT is applied to get an 8x8 block f of 
DCT coefficients. Each coefficient represents the 
amount of a particular spatial-frequency content of f. 5 
The lowest frequency coefficient (also called the DC 
coefficient) is f(0.0) while the highest frequency coef- 
ficient is f(7,7). 

2. An 8x8 matrix of integers Q. called the quantization 
table, is used to quantized the coefficients in f to form 20 
f c . For notational convenience, we number the 64 
entries in each 8x8 image block and each 8x8 block of 
DCT coefficients in raster-scan order, and use this 
ordering to refer to individual entries in the various 
blocks. Thus. f(u,v) is referred to as ff8u+v]. Using this 25 
notation, 

/alnMnl/ZClnl 0SnS63 


i 63 

We) = -&- 

The distortion D(I,Q) between an image I and its approxi- 
mation T due to quantization by Q is similarly defined as the 
mean-squared-error between corresponding pixel values. 
Clearly. D(I.Q) is the mean value of all D(f.Q) over all 
constituent blocks f in I. The distortion D(I,Q) is used to 
judge the “quality" of T. in quality measures such as 
Signal-to-Noise Ratio (SNR) and Peak Signal-to-Noise 
Ratio (PSNR): 

i/2 

PSNR = 10 log 10 ft ; g| 

SNR = 10 log 10 ftlo )' 

Here S is the mean-squared pixel value over I. Higher 
distortion implies poorer quality and vice versa. As entries 
in Q are made higher. D(I.Q) tends to increase. 

The distortion between the DCT coefficient-block f and its 
approximation f resulting from quantization by Q. is 

1 63 - 

D<f,Q) = -z r ^(Anl-AUf 


where. // represents division followed by rounding to 30 ^ designing rat e-distortion-optimal quantization tables, the 
the nearest integer. For example, for any following property of the DCT can be exploited: 


b > Q<a/A? = 


|_f +a5 J- 


ifakO 

ifa<0 


It will be appreciated that LxJ calculates the largest 
integer^ x. 

3. The block f G is entropy-coded (e.g., with Huffman 
codes) to exploit similarities across blocks to give the 
compressed block E(f e ). The sequence of these com- 
pressed blocks forms the compressed image. 

The decompression process reverses these steps as fol- 
lows: 

1. Each entropy -coded block E(f e ) is decoded to get the 
corresponding block of quantized coefficients 

2. Dequantization is done to construct the block f . as 
follows: 

/|nl=/aln|'filnl 0SnS63 


35 OOIGHtfe ) 


That is. the mean squared error in the pixel-domain is the 
same as the mean squared error in DCT-domain. This can be 
seen as follows: 


40 


45 


DCT(f)=f 


IDCT(f) =/ 

DCM=f 

Dcnf-f) -f-f 

63 63 - 

£««]-/[ n\f= I W-W 


using linearity of DCT 
since DCT preserves L 2 norms 


D(f,Q) = D{f,Q) 


50 This implies that D(I.Q) can be split as a sum into distortions 
in various DCT coefficients. Let D„(I. q) be defined as: 


3. The two-dimensional Inverse Discrete Cosine Trans- 
form (IDCT). is applied f to get the decompressed 
image block f. These decompressed blocks form the 
decompressed image T. 

The lossiness of the compression is essentially because of 
the quantization step (f— »f G ), as in general. 

Q\n\)Q\n\*j\n\ 


= Mean( (/)n| 

Where f|n)=(f|n]//q)-q and the mean is taken over all the 
blocks in the image. Then 

63 

D(1,Q)= I ZW.GM) 


55 


60 


This causes differences in pixel values between the original 
image block f and its approximation, the decompressed 
image block f . The mean-squared-error between f and f is 65 
known as distortion in f caused by Q and is denoted as 
D(f.Q). 


Since the focus of this discussion is on the effects of 
different quantization tables on a given image, the distortion 
D(I,Q) can be denoted simply as D(Q). Similarly, the 
distortion D„( I.q) in the n' h coefficient can be denoted by 
D„(q). Then. 
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0(2)= I AtGM) 

n=0 

This decomposition of D(Q) into contributions from indi- 5 
vidual coefficients is important as it simplifies the task of 
minimizing D(Q) to that of minimizing a sum, each of 
whose components depends on just one entry in Q. The rate 
(or compressed size) resulting from Q can also be split 
similarly into a sum. These two properties of the DCT allow 10 
the problem of optimizing the rate-distortion tradeoff to be 
solved using a dynamic programming approach for mini- 
mizing two sums. 

The degree of compression achieved is usually expressed 
in terms of the rate of the compressed image, which is the 15 
number of bits used per pixel: 

size of compressed image in bits 
~~ number of pixels in the image 

20 

Low rates are achieved when the quantized blocks f Q have 
similar entries (low entropy). The most common case is that 
of a coefficient being quantized to zero. The more zeros there 
are in f„ the fewer bits it would take to store it Thus, 

w 25 

increasing the entries of Q tends to decrease the rate. We 
denote the rate resulting with the use of a particular quan- 
tization table, Q, as R(Q). 

DCT has the nice property of being very close to the 
Karhunen-Loeve-Hotelling transform, a transform that pro- 
duces uncorrelated coefficients. The lack of correlation 
between coefficients allows the rate to be decomposed as a 
sum of contributions from individual coefficients. It will be 
appreciated by those skilled in the art that the coefficient- 
wise average of entropies of the quantized DCT coefficients 
is a very good estimate of the rate resulting from two-pass J 
Huffman coding of run lengths. This means that R(Q) can be 
approximated as a sum of rates of individuals coefficients. 

Let R„(q) be defined as 


6 

We call the quantization tables Q that satisfy these condi- 
tions (for some A or B) RD-optimal. 

RD-OPT takes an image I as input and produces 
RD-optimal quantization tables for a wide range of rates and 
distortions. The contributions to rate and distortion of indi- 
vidual coefficients, R„(Q[n]) and D„(Q[n]), respectively, just 
depend on the entry Q[n] of Q. RD-OFT first calculates 
(R„(Q[nJ) and D„(Q[ n]> for each possible value of Q|n], and 
then uses a dynamic programming approach to minimize 
sums of R„(Q[n]) against sums of D„(Q|n|). To calculate 
R„(q) and D„(q) for each possible q, a preliminary pass 
through the image is run to gather DCT statistics which are 
used in a novel way. 

RD-OPT, as shown in FIG. 2 steps 200 through 210, has 
several main steps. An image I of width W and height H. 
with pixel values in the range [0 ... M] is input. DCT 
statistics for the image are gathered 202. The statistics are 
used 204 to calculate R n (q) and D„(q) for each possible q. 
Dynamic programing is used to optimize 206 R(Q) against 
D(Q) and generate 208 RD-optimal DCT quantization tables 
Q based on an optimized R„(q). 

The tasks to the gathering DCT statistics step are further 
detailed in flowchart elements 212 through 222 shown in 
HG. 3. These DCT statistics for the image should answer the 
questions: 

1. How many times does the a’ h coefficient get quantized 
to value v when Q[n]=q? 

2. What is the mean-squared error for the n th coefficient 
when Q[n]=q? 

For any real number c, let 

{ |_2cJ ifc S 0 

-l-ad ifc< o 

It can be shown that for any integer q£l. 


Entropy) (fin ]//,})} 


40 


c//^=Bucket(cy/2^. 


Where the^ entropy is measured over all the blocks in the 
image. If (f[n]//q) takes the value v in a fraction p^>0 of all 
blocks f, then this entropy is -I v p,, log, p„. Then. 

63 

m>-~ XRJ.Q[n]). 


Thus, R(Q) can also be decomposed into a sum of contri- 
butions from individual coefficients, just like D(Q). 

It is desirable to have a low rate (high compression) and 
a low distortion (high quality). However, varying Q has 
opposite effects on D(Q) and R(Q). The distortion D(Q) 55 
tends to increase and the rate R(Q) tends to decrease as the 
entries in Q are made larger. The tradeoff between D(Q) and 
R(Q) is different to different images. The following discus- 
sion is directed to RD-OPT. an algorithm for constructing 
quantization tables with optimal rate-distortion tradeoffs to ^ 
a given image. 

The problem of choosing Q to optimize the rate-distortion 
tradeoff for a given image, L can be stated in two ways: 

1. Given a target distortion A, find Q such that D(Q)=A 

and the rate R(Q) is minimized. 65 

2. Given a target rate B bits per pixel (bpp), find Q such 
that R(Q)SB and the distortion D(Q) is minimized. 


Hence, it suffices to gather statistics by counting the number 
of times each DCT coefficient takes a value in a particular 
bucket, as this count can then be used to calculate the 
number of times a particular quantized coefficient takes a 
particular value. The unquantized coefficient value itself can 
be approximated to within +0.25. It will be appreciated by 
those skilled in the art that bucket (0) which corresponds to 
values in the interval (-0.5 ,0.5) must be split into two 
buckets, but this detail is ignored for the following discus- 
sion for added clarity and simplicity. This simplification 
introduces some error to the calculations, but the simplified 
calculations are still monotonically related to the more 
accurate calculations and as such will produce nearly the 
same result. 

Procedure GatherStats. as shown in Table 1 as a series of 
pseudo-software code, fills an array OccursCount|0 . . . 
63H-2VMAX . . . 2VMAX]. The constant VMAX is the 
maximum absolute value any DCT coefficient can take (i.e., 
for 1-byte samples. M=255 and VMAX=2048). The value 
OccursCount[n][v] at the end is the number of blocks where 
the n th DCT coefficient c n is such that Bucket(c„)=v. 

GatherStats works by going through each 8x8 block f in 
1 214 and calculating its Discrete Cosine Transform g 216. 
For each coefficient g[n|. the count OccursCount[n]lBucket 
(gin]) | is incremented by one 218. 
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TABLE 1 


Procedure GathersStats 
Input: Image I 

Output: Array OccursCountfO ... 63] [-2VMAX . . . 2VMAX] 

1 . Initialize Occur sCount to 0 eveiywhere 

2. For each 8x8 block f in I 

3. g := DCT(f) 

4. For n := 0 to 63 

5. v := Bucket(g[n}) 

6. OccursCoimt|n] [v] -H- 


The step of calculating R„(q) and D„(q) is further detailed 
in FIGS. 4 and 5 in steps 224 through 284. Let the possible 
range of values of any quantization table entry q be 
l^qgMAXQ. Let R[0 . . . 63][1 . . . MAXQ] and 
D[0 . . . 631[1 . . . MAXQj be two-dimensional arrays. The 
task of this step is to fill these arrays using the array 
OccursCount such that. 

R[n\[q\=R„(q) 

D\n\[q^DJq) 

This is accomplished by procedures FillR and FillD. 

FillR, shown in FIG. 4 and detailed in pseudo- software 
code in Table 2. fills the array R by calculating the entropy 
242 of the n"* coefficient when quantized by q, for all n and 
q. For each possible quantized value (QuantizedVal). the 
variable Count is used to compute the number of times the 
n" 1 coefficient gets quantized (by q) to QuantizedVal 232. 
Count is simply the sum of all OccursCount |n||v] such that 
v//(2q) is equal to QuantizedVal. If F is the total number of 
blocks in I. then, the entropy is calculated as: 

TABLE 2 


Procedure FillR 

Input: Array OccursCount JO ... 63] [-2VMAX . . . 2VMAX] 
Output: Array R|0 . . . 63] [1 . . . MAXQ] 

1. F := Number of 8x8 blocks in the image 

2. For n := 0 to 63 

3. For q := 1 to MAXQ 

4. Entropy = 0 

5. For QuantizedVal := (-VMAX) // q to VMAX // q 

/* QuantizedVal is the quantized value */ 

6. Count := 0 1* Count is the # of times the 

value QuantizedVal occurs */ 

7. For each v such that v // (2q) = 

QuantizedVal 

8. Count := Count + OccursCountfn] [v] 

9. Probab := Count/F 

10. If (Probab > 0) then 

11. Entropy := Entropy - (Probab * 

Log2 (Probab)) 

12. R In ] [q] := Entropy/64.0 


FillD. shown in FIG. S and detailed in pseudo- software 
code in Table 3. fills the array D by calculating the error in 
quantizing the n ,h coefficient by q. for each n and q. For each 
integer v in the range -2VMAX . . . 2VMAX, the n' h 
coefficient gets quantized to the value QuantizedVal (=v// 
(2q)) in OccursCount|n][v| blocks. The actual (unquantized) 
value of the coefficient in each of these blocks is estimated 
by the variable OrginalVal 262. to within ±0.25. Hence for 
each v. the error is incremented by OccursCount] n][v| times 
(OriginalVal-q times QuantizedVal) 2 262. 


TABLE 3 


Procedure FillD 

Input: Amy OccursCountfO ... 63] (-2VMAX . . . 2 VMAX] 
Output: Array D(0 . . . 63] [1 . . . MAXQ] 

1. N := Number of pixels in the image 

2. For n := 0 to 63 

3. For q := q to MAXQ 

4. Dfn] [q] := 0 

5. For v := -2VMAX to 2 VMAX 

*/ OriginalVal is the original coefficient 
value, within 0.25 */ 

6. OriginalVal = v/2.0 + ((v < 0) ? -0.25 : 0.25) 

/* QuantizedVal is the quantized value */ 

7. QuantizedVal — y II (2g) 

8. Error := OccursCount [n] fvj * S quare( Original Val 

- q*QuantizedVal) 

9. Dfn] fqj ;= Dfn] [q] + Eitot 
10. Dfn] [q] := Dfn] fqj/N 


The step of finding RD-optimal Q is further detailed in 
FIG. 6 in steps 286 through 310. The arrays R[ J and 
D[. . . ]|. . . ] are used in this step to find rate-distortion- 
optimal quantization tables. A dynamic programing (DP) 
approach is used. For this approach either R(Q) or D(Q) 
must be discretized to integral values so that it can be used 
as an index in the DP table. It will be appreciated by those 
skilled in the art that this will cause errors in the quantity 
discretized. 

Let BPPSCALE be a large integer constant. R(Q) is 
discretized to be integral by multiplying each R|n][q] with 
BPPSCALE and rounding off. This can be done in Proce- 
dure FillR itself, by inserting Ihe following pseudo- software 
code lines after line 12 as detailed in Table 4. Alternatively, 
this can be done by the same pseudo software code lines by 
inserting them before line 1 in the Procedure FillLeastD. 

TABLE 4 


For n := 0 to 63 

For q := 1 to MAXQ 

R|n] [q] := (R[n| |q] ‘BPPSCALE) III 


For the rest of this section, only discretized values will be 
referred to when referring to R(Q). and R„(q)=R(n][q]. Let 
MAXRATE be the discretized value of the highest rate for 
which we are interested in finding an RD-optimal quantiza- 
tion table. Let LeastDjO . . . 63 1[ 0 . . . MAXRATE] be an 
array whose entries have the following definition: LeastD 
]n][sj is the least total distortion for coefficients numbered 0 
through n such that the total rate (for these coefficients) is 
exactly s. That is. LeastD|nj[sJ is the least value (over all Q) 
of 

2 


subject to the constraint 


The procedure starts with each entry in LeastD set to°c288 
and then fills the rows one by one. The key idea is described 
in Theorem 1. 
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Theorem 1: 


For each n, l n = 63, and each s, 0 £ s £ MAXRATE, 
let X(n, s) be the set 




1 S q £ MAXQ 

X(n. j) = 

D[n][ 9 | + Least£>ln-llPi 

s' = s- R[rt][q] 

j go 


10 

the 63rd row. To recover the desired quantization table Q 
from that point, procedure RecoverQ. shown below in Table 
6 in pseudo-software code and which corresponds to step 
308 in FIG. 6, is used. This procedure recovers the quanti- 
zation table Q for target rate s by setting Q[63] to QChoice 
[63][s] and then working its way up the rows as follows. For 
going from row number n to row number n-1. s is decre- 
mented by R|n)[Q [nj] and then Q| n— 1 ] is set to QChoice 
[n — 1 J[sJ. 


Then, 


TABLE 6 


( min X( n, s) 


if Xt'n, j) is non-empty 
otherwise 


Proof: Suppose X (n.s) is empty or min X (n.s) is «■. Then. 15 
for every q(l<q<MAXQ). either R(n][q|>s, or LeastD[n-l] 
[s-R[n][q]]=°°. In either case, the rate s cannot be achieved 
from coefficients 0 through n. implying LeastD[n][s]=°°. 

Now assume X(n.s) is non-empty and that d is the 
minimum value in X (n,s). achieved by setting Q[n] to q. 20 
Assume LeastD[n][s]=d’>d. Then the distortion d' must be 
achieved with some value say q 1 for Q[n], Let d"=d'-D[n] 

[q’]. Then the distortion d” must be achieved from coeffi- 
cients 0-through n-1, with their rate being exactly equal to 
s-R[n][q’]. But then. d"=LeastD[n-l][s-R[n][q'J], as other- 
wise d’ can be improved, contradicting d’=LeastD[n][s]. 25 
Hence, d'=D[n][q’]+LeastD[n-l][s-R[n][q'] implying d’eX 
(n.s). Thus d^d'. which contradicts d'>d. 

To recover the RD-optimal quantization table for any 
desired distortion or rate, another data structure 
QChoice [0 . . . 63] [0 . . . MAXRATE] is maintained. 30 
QChoicefn][s] stores the value q that gave the entry in 
LeastD[n][s]. 

The procedure FillLeastD. which fills the array LeastD. 
using Theorem 1 is shown below in Table 5 in pseudo- 
software code and corresponds to steps 288 to 306 in FIG. 35 
6 . Row number n is filled using row number n-1, D[n][. . . 

], and R[n][. . . ]. For each q and each s', D[n][q]+LeastD 
fn-l][s'] is compared with LeastD[n][s'+R[n][q]]. If the 
former is lesser, then it replaces the latter, with QChoice[n] 
|s'+R[n][q]| being set to q. ^ 

TABLE 5 


Procedure FillLeastD 

Input: Arrays D[0 . . . 63 ] [1 - . . MAXQ], 

R[0 . . . 63] [0 . . . MAXSIZE] 45 

Output: Arrays LeastD[0 . . . 63 ] [0 . . . MAXSIZE, 

Qcboice[0 . . . 63] [0 . . . MAXSIZE] 

/* Initializations */ 

1 . For n := 0 to 63 

2. For s := 0 to MAXSIZE 

3. LeastD{n] [s| := INFINITY 
/* Fill row number zero */ 

4. For q := 1 to MAXQ 

5. If (D[0] [q] < LeastD|0] [R[0] [q]] ) then 

6. LeastDfO] |R|0) JqlJ := D[0][q| 

7. QChoke|0| |R[0[ [q]| := q 
/* Main loop */ 

8. For n := 1 to 63 $5 

9. For q := 1 to MAXQ 

10. For s' := 0 to MAXRATE 

11. If (D[n| [q] + LeastDfn-l] [s'] <LeastD[n] [s' 

+ R|n| |q][) Then 

12. s := s' + R[n) [q| 

13. LeastD[n| [s] := D[n] (qj + LeastD[n-l] [s'] 60 

14. QChoice[n| [s] := q 


Now, if a total distortion requirement A is to be met, it is 
straightforward to find the least s such that LeastD[63][s] 
§A. Similarly, if a rate requirement B is to be met. it is easy 65 
to find s such that sSB and LeastD[63][s] is the minimum 
over all such s. Thus, in both cases, a starting point s is in 


Procedure RecoverQ 

Input: Array QChoice[0 . . . 63] |0 . . . MAXSIZE1; Target rate s 
Output: Quantization table Q[0 ... 63 1 

1. For n := 63 down to 0 

2. Q[n] := Qchoicefn] [s] 

3. s := s - R[n| [Q[n[J 


It will be appreciated by those skilled in the art that the 
GatherStats procedure should execute in about the time that 
is required to apply the DCT once to each block in the 
image. Also, the FillD and FillR procedures each run in less 
time than a constant times 64xMAXQxVMAX. 

In addition, the FillLeastD procedure runs in less time 
than a constant times 64xMAXQxMAXRATE. In any prac- 
tical implementation, this can be substantially reduced The 
loop range in line 9 of Table 5 can be made 1 to MAXQ(n) 
where MAXQ(n) is the minimum of MAXQ and the least 
value of Q[ n] that will make the n' h coefficient zero every- 
where in the image. The loop range in line 10 for s' can be 
made 1 to the last entry in the (n-1)'* row which is not 
marked °°. Further, if only one given target distortion or 
target rate is to be met. then the loops can be pruned to 
exclude cases which will clearly be outside the given 
specifications. Also, only two rows of the table LeastD need 
to be maintained at any point: the current row and the 
previous row. 

Further, the key idea is entirely symmetrical, in the sense 
that an array LeastR can be maintained instead of LeastD. As 
a result, a procedure FillLeastR analogous to FillLeastD 
would be used. In this case, distortion would need to be 
discretized to integral values. The choice of the procedure to 
be used depends on various factors such as the range of 
values spanned by distortion and rate, and the errors incurred 
by discretizing them to integers. 

The RD-OFT algorithm has another interesting feature in 
that the final results are independent of the order in which 
the coefficients are considered. This implies that it can be 
readily parallelized. The 64 rows can be pairwise combined, 
then the 32 “composite” rows can be pairwise combined, 
and so on. 

Another variation is to use weighted mean squared error 
instead of mean squared error by assigning different weights 
to errors in different frequencies in FillD. This might be used 
to give distortion in lower frequencies more importance. 

For better visual quality, it may be useful to do adaptive 
quantization which gives more bits for encoding regions in 
the image that are perceptually more significant. This is done 
in MPEG by scaling the quantization table up or down on a 
per- macroblock basis. Thus, for any block f. the quantization 
table used is Q times qscale^, where Q is a nominal quan- 
tization table and qscale^ is a factor that depends upon the 
macroblock containing f. The value of q scale,, should be 
chosen based upon characteristics such as texture, total 
energy, presence of edges, etc. However, qscaky does not 
depend upon Q. Hence, while gathering statistics (i.e.. 
procedure GatherStats) qscaie / can be determined for each 
block. The entry OccursCount[n][v] can be filled by setting 
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v to be the actual value of the n‘ h coefficient divided by 
qscalef for the block under consideration. Then. FillLeastD 
will optimize Q to give the best rate-distortion tradeoff for 
the adaptive quantization scheme. 

The present invention can be summarized in reference to 
FIG. 1 which is a block diagram of the preferred embodi- 
ment discrete cosine transform-based digital image com- 
pression and decompression system 100. The system 100 
includes a transmitting device which receives input digital 
image data 102. generates compressed digital image data 
116. and transmits, through a transmitter 112, the com- 
pressed digital image data 114 and a rate-distortion-optimal 
quantization table 114 over a communication channel 118. 
The transmitting device generates compressed digital image 
data 116 by utilizing a discrete cosine transformer 104 which 
divides the input digital image data 102 into fixed length 
image blocks and generates a corresponding block of dis- 
crete cosine transform coefficients for each fixed length 
image block by applying a discrete cosine transform to each 
fixed length image block. 

A quantization table generator 106 (i.e., a digital image 
compression preprocessor 106) is operatively coupled to the 
discrete cosine transformer 104 to generate the rate- 
distortion-optimal quantization table by determining dis- 
crete cosine transform statistics from each block of discrete 
cosine transform coefficients, calculating an image distor- 
tion array and a rate of image compression array based upon 
the discrete cosine transform statistics for each possible 
quantization value, and optimizing the rate of image com- 
pression array against the image distortion array to derive 
the rate-distortion-optimal quantization table. 

A quantizer 108 is operatively coupled to the discrete 
cosine transformer 104 and the quantization table generator 
106 to generate a corresponding block of quantized discrete 
cosine transform coefficients for each block of discrete 
cosine transform coefficients by quantizing each block of 
discrete cosine transform coefficients with the rate- 
distortion-optimal quantization table. An entropy coder 110 
is operatively coupled to the quantizer 108 to entropy 
encode (such as with Huffman codes which exploit simi- 
larities across the quantized discrete cosine transform 
coefficients) each block of quantized discrete cosine trans- 
form coefficients to form a corresponding compressed block 
such that a sequence of all of the compressed blocks forms 
the compressed digital image data 

It will be appreciated by those skilled in the art that the 
communication channel can take many forms including, but 
not limited to. a serial communication line, a parallel com- 
munication line, a data network, a public switched telephone 
network, a optically-encoded disk, a magnetically-encoded 
disk, and/or a radio communication link. 

A receiving device which receives, through a receiver 
124. the compressed digital image data 122 and the rate- 
distortion-optimal quantization table 122 from the commu- 
nication channel 118 and generates 120 decompressed digi- 
tal image data 132 is also a part of the system 100. The 
receiving device generates 120 decompressed digital image 
data 132 by utilizing a decoder 126 to decode each entropy- 
coded compressed block in the received compressed digital 
image data 122 into a corresponding block of quantized 
discrete cosine transform coefficients. A dequantizer 128 is 
operatively coupled to the decoder 126 to generate a corre- 
sponding block of discrete cosine transform coefficients for 
each block of quantized discrete cosine transform coeffi- 
cients based on the received rate-distortion-optirnal quanti- 
zation table 122. An inverse discrete cosine transformer 130 
is operatively coupled to the dequantizer 128 to generating 
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a corresponding decompressed block of fixed length digital 
image data for each block of discrete cosine transform 
coefficients by applying an inverse discrete cosine transform 
to each block of discrete cosine transform coefficients such 
5 that a sequence of all of the decompressed blocks forms the 
decompressed digital image data 132. 

The system 100 may also include an input device which 
is operatively coupled to the transmitting device for receiv- 
ing the input digital image data from a device external to the 
to digital image compression and decompression system 100. 
The input device may be a serial communications port, a 
parallel communications port, a data network interface (i.e.. 
to a computer data bus, ethernet, fiber-optic network or the 
like), a public switched telephone network (PSTN) adapter, 
15 a optically-encoded disk drive, a magnetically-encoded disk 
drive (e.g.. a floppy or hard disk drive), and/or a radio 
communication unit for a trunked radio network, cellular 
radio voice or data packet network, personal communication 
network, and/or satellite system. The transmitter 112 in 
20 system 100 may also transmit the compressed digital image 
data and the rate-distortion-optimal quantization table 
through the same sort of devices and the input devices. 

The quantization table generator 106 of system 100 
preferably is configured to generate the rate-distortion- 
25 optimal quantization table according to a particular algo- 
rithm. This algorithm includes determining discrete cosine 
transform statistics from input digital image data by dividing 
the input digital image data into fixed length image blocks, 
applying a discrete cosine transform to each fixed length 
30 image block to generate a corresponding block of discrete 
cosine transform coefficients for each fixed length image 
block, and counting a number of times each discrete cosine 
transform coefficient has a non-zero value for all of the fixed 
length image blocks divided out of the input digital image 
35 data. In addition, the algorithm includes calculating the rate 
of image compression array by calculating entropy of each 
discrete cosine transform coefficient when quantized by each 
possible quantization value as a function of the counted 
number of times each discrete cosine transform coefficient 
40 has a non-zero value. Also, the algorithm includes calculat- 
ing the image distortion array by calculating error produced 
when quantizing each discrete cosine transform coefficient 
by each possible quantization value as a function of the 
counted number of times each discrete cosine transform 
45 coefficient has a non-zero value. Further, the algorithm 
includes deriving the rate-distortion-optimal quantization 
table by calculating a least distortion array containing the 
least distortion produced for each discrete cosine transform 
coefficient at each rate of image compression as a function 
50 of the image distortion array and the rate of image com- 
pression array, calculating an array of quantization values 
which produced each entry in the least distortion array, and 
determining the rate-distortion-optimal quantization table at 
least as a function of the array of quantization values which 
55 produced each entry in the least distortion array. 

The quantization table generator 106 of system 100 may 
also include a mechanism for discretizing values in one of 
the image distortion array and the rate of image compression 
array to integral values prior to deriving the rate-distortion- 
60 optimal quantization table. 

The quantization table generator 106 of system 100 
preferably includes a mechanism for optimizing the rate of 
image compression array against the image distortion array 
with respect to a total image distortion requirement to derive 
65 the rate -distortion-optimal quantization table having a maxi- 
mal rate of image compression with a total image distortion 
less than the total image distortion requirement. 
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Alternatively, the quantization table generator 106 of 
system 100 may include a mechanism for optimizing the rate 
of image compression array against the image distortion 
array with respect to a rate of image compression require- 
ment to derive the rate-distortion-optimal quantization table 
having a minimal total image distortion with a rate of image 
compression greater than the rate of image compression 
requirement. 

Also, the quantization table generator 106 of system 100 
preferably includes apparati for configuring the rate- 
distortion-optimal quantization table for use with various 
image compression standard algorithms such that the digital 
image compression and decompression system 100 utilizes 
any one of the following image compression standard 
algorithms, including: a joint photographic experts group 
(JPEG) image compression standard algorithm, a motion 
picture experts group (MPEG) image compression standard 
algorithm, and a H.261 video compression standard algo- 
rithm 

It will be appreciated by those skilled in the art that the 
system 100 can be subdivided into distinct parts and imple- 
mented separately based on function. For example, a single 
discrete cosine transform-based digital image compression 
device 116 could be produced. Similarly, the functions of the 
quantization table generator 106 may be implemented as a 
digital image compression preprocessor 106. 

Also, the system 100 and its parts may be described in 
reference to methods for: generating a rate -distortion- 
optimal quantization table, using discrete cosine transform- 
based digital image compression, and operating a discrete 
cosine transform-based digital image compression and 
decompression system 

Although the invention has been described and illustrated 
with a certain degree of particularity, it is understood that the 
present disclosure of embodiments has been made by way of 
example only and that numerous changes in the arrangement 
and combination of parts as well as steps may be reseated to 
by those skilled in the art without departing from the spirit 
and scope of the invention as claimed. For example, pseudo- 
software code described above is provided for discussion 
purposes only. This underlying algorithm could be imple- 
mented in other pseudo-software code without departing 
from the scope and spirit of the present invention as claimed. 
In addition, the pseudo-software code could be implemented 
in any computer software language or set in firmware 
without departing from the spirit and scope of the invention 
as claimed. 

What is claimed is: 

1. A digital image compression preprocessor for use in a 
discrete cosine transform-based digital image compression 
device, the preprocessor comprising: 

(a) gathering means for determining discrete cosine trans- 
form statistics from input digital image data, wherein 
the discrete cosine transform statistics include counts 
related to the number of times each of a plurality of 
discrete cosine transform coefficients derived from the 
input digital image data takes a value; 

(b) computing means, operatively coupled to the gather- 
ing means, for calculating an image distortion array and 
a rate of image compression array based upon the 
discrete cosine transform statistics for each of a plu- 
rality of possible quantization values; and 

(c) quantization table generating means, operatively 
coupled to the computing means, for optimizing the 
rate of image compression array against the image 
distortion array to derive a rate-distortion-optimal 
quantization table. 
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2. The digital image compression preprocessor of claim 1 
wherein the gathering means comprises means for determin- 
ing discrete cosine transform statistics from input digital 
image data by dividing the input digital image data into fixed 

5 length image blocks, applying a discrete cosine transform to 
each fixed length image block to generate a corresponding 
block of discrete cosine transform coefficients for each fixed 
length image block, and counting a number of times each 
discrete cosine transform coefficient has a non-zero value for 
to all of the fixed length image blocks divided out of the input 
digital image data. 

3. The digital image compression preprocessor of claim 2 
wherein the computing means comprises means for calcu- 
lating the rate of image compression array by calculating 

is entropy of each discrete cosine transform coefficient when 
quantized by each possible quantization value as a function 
of the counted number of times each discrete cosine trans- 
form coefficient has a non-zero value. 

4. The digital image compression preprocessor of claim 3 
20 wherein the computing means comprises means for calcu- 
lating the image distortion array by calculating error pro- 
duced when quantizing each discrete cosine transform coef- 
ficient by each possible quantization value as a function of 
the counted number of times each discrete cosine transform 

25 coefficient has a non-zero value. 

5. The digital image compression preprocessor of claim 4 
wherein the quantization table generating means comprises 
means for deriving the rate-distortion-optimal quantization 
table by calculating a least distortion array containing the 

30 least distortion produced for each discrete cosine transform 
coefficient at each rate of image compression as a function 
of the image distortion array and the rate of image com- 
pression array, calculating an array of quantization values 
which produced each entry in the least distortion array, and 
35 determining the rate-distortion-optimal quantization table at 
least as a function of the array of quantization values which 
produced each entry in the least distortion array. 

6. The digital image compression preprocessor of claim 5 
wherein the quantization table generating means comprises 

40 means for discretizing values in one of the image distortion 
array and the rate of image compression array to integral 
values prior to the deriving of the rate-distortion-optimal 
quantization table. 

7. The digital image compression preprocessor of claim 1 
45 wherein the quantization table generating means comprises 

means for optimizing the rate of image compression array 
against the image distortion array with respect to a total 
image distortion requirement to derive a rate-distortion- 
optimal quantization table having a maximal rate of image 
50 compression with a total image distortion less than the total 
image distortion requirement. 

8. The digital image compression preprocessor of claim 1 
wherein the quantization table generating means comprises 
means for optimizing the rate of image compression array 

55 against the image distortion array with respect to a rate of 
image compression requirement to derive a rate-distortion- 
optimal quantization table having a minimal total image 
distortion with a rate of image compression greater than the 
rate of image compression requirement. 

60 9. The digital image compression preprocessor of claim 1 

wherein the digital image compression device utilizes joint 
photographic experts group (JPEG) image compression 
standard algorithms and wherein the quantization table 
generating means comprises means for configuring the rate- 
65 distortion-optimal quantization table for use with the joint 
photographic experts group (JPEG) image compression 
standard algorithms. 
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10. The digital image compression preprocessor of claim 
1 wherein the digital image compression device utilizes 
motion picture experts group (MPEG) image compression 
standard algorithms and wherein the quantization table 
generating means comprises means for configuring the rate- 
distortion-optimal quantization table for use with the motion 
picture experts group (MPEG) image compression standard 
algorithms. 

11. The digital image compression preprocessor of claim 
1 wherein the digital image compression device utilizes 
H.261 video compression standard algorithms and wherein 
the quantization table generating means comprises means 
for configuring the rate-distortion-optimal quantization table 
for use with the H.261 video compression standard algo- 
rithms. 

12. A discrete cosine transform-based digital image com- 
pression device, comprising: 

(a) discrete cosine transformer means for dividing input 
digital image data into fixed length image blocks and 
generating a corresponding block of discrete cosine 
transform coefficients for each fixed length image block 
by applying a discrete cosine transform to each fixed 
length image block; 

(b) quantization table generator, operatively coupled to 
the discrete cosine transformer means, which generates 
a rate-distortion-optimal quantization table, compris- 
ing: 

(i) gathering means for determining discrete cosine 
transform statistics from each block of discrete 
cosine transform coefficients, wherein the discrete 
cosine transform statistics include counts related to 
the number of times each discrete cosine transform 
coefficient from the blocks takes a value; 

(ii) computing means, operatively coupled to the gath- 
ering means, for calculating an image distortion 
array and a rate of image compression array based 
upon the discrete cosine transform statistics for each 
of a plurality of possible quantization values; and 

(iii) quantization table generating means, operatively 
coupled to the computing means, for optimizing the 
rate of image compression array against the image 
distortion array to derive the rate-distortion-optimal 
quantization table; 

(c) quantizing means, operatively coupled to the discrete 
cosine transformer means and the quantization table 
generating means, for generating a corresponding block 
of quantized discrete cosine transform coefficients for 
each block of discrete cosine transform coefficients by 
quantizing each block of discrete cosine transform 
coefficients with the rate-distortion-optimal quantiza- 
tion table; and 

(d) entropy coding means, operatively coupled to the 
quantizing means, for entropy coding each block of 
quantized discrete cosine transform coefficients to form 
a corresponding compressed block such that a sequence 
of all of the compressed blocks forms the compressed 
digital image data. 

13. The digital image compression device of claim 12 
wherein the entropy coding means comprises means for 
entropy coding each block of quantized discrete cosine 
transform coefficients with Huffman codes to exploit simi- 
larities across the quantized discrete cosine transform coef- 
ficients. 

14. The digital image compression device of claim 12 
further comprising input means, operatively coupled to the 
discrete cosine transformer means, for receiving the input 
digital image data from a device external to the digital image 
compression device. 
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15. The digital image compression device of claim 14 
wherein the device external to the digital image compression 
device is selected from the group consisting of a serial 
communications port, a parallel communications port, a data 

5 network interface, a public switched telephone network 
adapter, an optically-encoded disk drive, a magnetically- 
encoded disk drive, and a radio communication unit. 

16. The digital image compression device of claim 12 
further comprising transmitting means, operatively coupled 
to the entropy coding means, for transmitting the com- 
pressed digital image data and the rate-distortion-optimal 
quantization table to a device external to the digital image 
compression device. 

17. The digital image compression device of claim 16 
wherein the device external to the digital image compression 

15 device is selected from the group consisting of a serial 
communications port, a parallel communications port, a data 
network interface, a public switched telephone network 
adapter, an optically-encoded disk drive, a magnetically- 
encoded disk drive, and a radio communication unit. 

20 18. The digital image compression device of claim 12 

wherein; 

(a) the quantization table generator gathering means com- 
prises means for determining discrete cosine transform 
statistics from input digital image data by dividing the 

25 input digital image data into fixed length image blocks, 
applying a discrete cosine transform to each fixed 
length image block to generate a corresponding block 
of discrete cosine transform coefficients for each fixed 
length image block, and counting a number of times 

30 each discrete cosine transform coefficient has a non- 
zero value for all of the fixed length image blocks 
divided out of the input digital image data; 

(b) the quantization table generator computing means 
comprises means for calculating the rate of image 

35 compression array by calculating entropy of each dis- 
crete cosine transform coefficient when quantized by 
each possible quantization value as a function of the 
counted number of times each discrete cosine transform 
coefficient has a non-zero value; 

40 (c) the quantization table generator computing means 

comprises means for calculating the image distortion 
array by calculating error produced when quantizing 
each discrete cosine transform coefficient by each pos- 
sible quantization value as a function of the counted 

45 number of times each discrete cosine transform coef- 
ficient has a non-zero value; and 

(d) the quantization table generator quantization table 
generating means comprises means for deriving the 
rate-distortion-optimal quantization table by calculat- 

50 ing a least distortion array containing the least distor- 
tion produced for each discrete cosine transform coef- 
ficient at each rate of image compression as a function 
of the image distortion array and the rate of image 
compression array, calculating an array of quantization 

55 values which produced each entry in the least distortion 
array, and determining the rate-distortion-optimal 
quantization table at least as a function of the array of 
quantization values which produced each entry in the 
least distortion array. 

60 19. The digital image compression device of claim 18 

wherein the quantization table generator quantization table 
generating means comprises means for discretizing values in 
one of the image distortion array and the rate of image 
compression array to integral values prior to deriving the 

65 rate-distortion-optimal quantization table. 

20. The digital image compression device of claim 12 
wherein the quantization table generator quantization table 
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generating means comprises means for optimizing the rate 
of image compression array against the image distortion 
array with respect to a total image distortion requirement to 
derive a rate-distortion-optimal quantization table having a 
maximal rate of image compression with a total image 
distortion less than the total image distortion requirement. 

21. The digital image compression device of claim 12 
wherein the quantization table generator quantization table 
generating means comprises means for optimizing the rate 
of image compression array against the image distortion 
array with respect to a rate of image compression require- 
ment to derive a rate-distortion-optimal quantization table 
having a minimal total image distortion with a rate of image 
compression greater than the rate of image compression 
requirement. 

22. The digital image compression device of claim 12 
wherein the quantization table generator quantization table 
generating means comprises means for configuring the rate- 
distortion-optimal quantization table for use with image 
compression standard algorithms such that the digital image 
compression device utilizes an image compression standard 
algorithm selected from the group consisting of a joint 
photographic experts group (JPEG) image compression 
standard algorithm, a motion picture experts group (MPEG) 
image compression standard algorithm, and a H.261 video 
compression standard algorithm. 

23. A discrete cosine transform-based digital image com- 
pression and decompression system, comprising: 

(a) a transmitting device which receives input digital 
image data, generates compressed digital image data, 
and transmits the compressed digital image data and a 
rate-distortion-optimal quantization table over a com- 
munication channel, the transmitting device compris- 
ing: 

(i) discrete cosine transformer means for dividing the 
input digital image data into fixed length image 
blocks and generating a corresponding block of 
discrete cosine transform coefficients for each fixed 
length image block by applying a discrete cosine 
transform to each fixed length image block; 

(ii) quantization table generating means, operatively 
coupled to the discrete cosine transformer means, for 
generating the rate-distortion-optimal quantization 
table by determining discrete cosine transform sta- 
tistics from each block of discrete cosine transform 
coefficients, wherein the discrete cosine transform 
statistics include counts related to the number of 
times each discrete cosine transform coefficient from 
the blocks takes a value, calculating an image dis- 
tortion array and a rate of image compression array 
based upon the discrete cosine transform statistics 
for each of a plurality of possible quantization 
values, and optimizing the rate of image compres- 
sion array against the image distortion array to derive 
the rate-distortion-optimal quantization table; 

(iii) quantization means, operatively coupled to the 
discrete cosine transformer means and the quantiza- 
tion table generating means, for generating a corre- 
sponding block of quantized discrete cosine trans- 
form coefficients for each block of discrete cosine 
transform coefficients by quantizing each block of 
discrete cosine transform coefficients with the rate- 
distortion-optimal quantization table; and 

(iv) entropy coding means, operatively coupled to the 
quantizing means, for entropy coding each block of 
quantized discrete cosine transform coefficients to 
form a corresponding compressed block such that a 
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sequence of all of the compressed blocks forms the 
compressed digital image data; and 

(b) a receiving device which receives the compressed 
digital image data and the rate-distortion-optimal quan- 
5 tization table from the communication channel and 
generates decompressed digital image data, the receiv- 
ing device comprising: 

(i) decoding means for decoding each entropy-coded 
compressed block in the received compressed digital 

10 image data into a corresponding block of quantized 

discrete cosine transform coefficients; 

(ii) dequantizing means, operatively coupled to the 
decoding means, for generating a corresponding 
block of discrete cosine transform coefficients for 

15 each block of quantized discrete cosine transform 

coefficients based on the received rate-distortion- 
optimal quantization table; and 

(iii) inverse discrete cosine transformer means, opera- 
tively coupled to the dequantizing means, for gen- 

20 crating a corresponding decompressed block of fixed 

length digital image data for each block of discrete 
cosine transform coefficients by applying an inverse 
discrete cosine transform to each block of discrete 
cosine transform coefficients such that a sequence of 
25 all of the decompressed blocks forms the decom- 

pressed digital image data. 

24. The digital image compression and decompression 
system of claim 23 wherein the communication channel is 
selected from the group consisting of a serial communica- 

30 tion line, a parallel communication line, a data network, a 
public switched telephone network, an optically-encoded 
disk, a magnetically-encoded disk, and a radio communica- 
tion link. 

25. The digital image compression and decompression 
35 system of claim 23 further comprising input means, opera- 
tively coupled to the transmitting device, for receiving the 
input digital image data from a device external to the digital 
image compression and decompression system. 

26. The digital image compression and decompression 
40 system of claim 25 wherein the input means is selected from 

the group consisting of a serial communications port, a 
parallel communications port, a data network interface, a 
public switched telephone network adapter, an optically- 
encoded disk drive, a magnetically-encoded disk drive, and 
45 a radio communication unit. 

27. The digital image compression and decompression 
system of claim 23 wherein the quantization table generating 
means is configured to generate the rate-distortion-optimal 
quantization table by: 

50 (a) determining discrete cosine transform statistics from 

input digital image data by dividing the input digital 
image data into fixed length image blocks, applying a 
discrete cosine transform to each fixed length image 
block to generate a corresponding block of discrete 
55 cosine transform coefficients for each fixed length 

image block, and counting a number of times each 
discrete cosine transform coefficient has a non-zero 
value for all of the fixed length image blocks divided 
out of the input digital image data; 

60 (b) calculating the rate of image compression array by 

calculating entropy of each discrete cosine transform 
coefficient when quantized by each possible quantiza- 
tion value as a function of the counted number of times 
each discrete cosine transform coefficient has a non- 
65 zero value; 

(c) calculating the image distortion array by calculating 
error produced when quantizing each discrete cosine 
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transform coefficient by each possible quantization 
value as a function of the counted number of times each 
discrete cosine transform coefficient has a non-zero 
value; and 

(d) deriving the rate-distorti on-optimal quantization table 
by calculating a least distortion array containing the 
least distortion produced for each discrete cosine trans- 
form coefficient at each rate of image compression as a 
function of the image distortion array and the rate of 
image compression array, calculating an array of quan- 
tization values which produced each entry in the least 
distortion array, and determining the rate-distortion- 
optimal quantization table at least as a function of the 
array of quantization values which produced each entry 
in the least distortion array. 

28. The digital image compression and decompression 
system of claim 27 wherein the quantization table generating 
means comprises means for discretizing values in one of the 
image distortion array and the rate of image compression 
array to integral values prior to deriving the rate-distortion- 
optimal quantization table. 

29. The digital image compression and decompression 
system of claim 23 wherein the quantization table generating 
means comprises means for opti mizin g the rate of image 
compression array against the image distortion array with 
respect to a total image distortion requirement to derive a 
rate -distortion-optimal quantization table having a maximal 
rate of image compression with a total image distortion less 
than the total image distortion requirement. 

30. The digital image compression and decompression 
system of claim 23 wherein the quantization table generating 
means comprises means for optimizing the rate of image 
compression array against the image distortion array with 
respect to a rate of image compression requirement to derive 
a rate -distortion-optimal quantization table having a mini- 
mal total image distortion with a rate of image compression 
greater than the rate of image compression requirement 

31. The digital image compression and decompression 
system of claim 23 wherein the quantization table generating 
means comprises means for configuring the rate-distortion- 
optimal quantization table for use with image compression 
standard algorithms such that the digital image compression 
and decompression system utilizes an image compression 
standard algorithm selected from the group consisting of a 
joint photographic experts group (JPEG) image compression 
standard algorithm, a motion picture experts group (MPEG) 
image compression standard algorithm, and a H.261 video 
compression standard algorithm. 

32. A method for generating a rate-distortion-optimal 
quantization table for use in a discrete cosine transform- 
based digital image compression process, the method com- 
prising the steps of: 

(a) determining discrete cosine transform statistics from 
input digital image data, wherein the discrete cosine 
transform statistics include counts related to the num- 
ber of times each of a plurality of discrete cosine 
transform coefficients derived from the input digital 
image data takes a value; 

(b) calculating an image distortion array and a rate of 
image compression array based upon the discrete 
cosine transform statistics for each of a plurality of 
possible quantization values; and 

(c) optimizing the rate of image compression array against 
the image distortion array to derive the rate-distortion- 
optimal quantization table. 

33. The method of claim 32 wherein the determining step 
comprises determining discrete cosine transform statistics 
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from input digital image data by dividing the input digital 
image data into fixed length image blocks, applying a 
discrete cosine transform to each fixed length image block to 
generate a corresponding block of discrete cosine transform 
5 coefficients for each fixed length image block, and counting 
a number of times each discrete cosine transform coefficient 
has a non-zero value for all of the fixed length image blocks 
divided out of the input digital image data. 

34. The method of claim 33 wherein the calculating step 
comprises calculating the rate of image compression array 
by calculating entropy of each discrete cosine transform 
coefficient when quantized by each possible quantization 
value as a function of the counted number of times each 
discrete cosine transform coefficient has a non-zero value. 

35. The method of claim 34 wherein the calculating step 
15 comprises calculating the image distortion array by calcu- 
lating eiror produced when quantizing each discrete cosine 
transform coefficient by each possible quantization value as 
a function of the counted number of times each discrete 
cosine transform coefficient has a non-zero value. 

20 36. The method of claim 35 wherein the optimizing step 

comprises deriving the rate-distortion-optimal quantization 
table by calculating a least distortion array containing the 
least distortion produced for each discrete cosine transform 
coefficient at each rate of image compression as a function 
of the image distortion array and the rate of image com- 
pression array, calculating an array of quantization values 
which produced each entry in the least distortion array, and 
determining the rate-distortion-optimal quantization table at 
least as a function of the array of quantization values which 
produced each entry in the least distortion array. 

37. The method of claim 35 wherein the optimizing step 
comprises discretizing values in one of the image distortion 
array and the rate of image compression array to integral 
values prior to the step of deriving the rate-distortion- 
optimal quantization table. 

38. The method of claim 32 wherein the optimizing step 
comprises optimizing the rate of image compression array 
against the image distortion array with respect to a total 
image distortion requirement to derive the rate-distortion- 
optimal quantization table having a maximal rate of image 
compression with a total image distortion less than the total 
image distortion requirement 

39. The method of claim 32 wherein the optimizing step 
comprises optimizing the rate of image compression array 
against the image distortion array with respect to a rate of 
image compression requirement to derive the rate- 
distortion-optimal quantization table having a minimal total 
image distortion with a rate of image compression greater 
than the rate of image compression requirement. 

40. The method of claim 32 wherein the method is 
modified such that the generated rate-distortion-optimal 
quantization table is configured for use with an image 
compression standard algorithm selected from the group 
consisting of a joint photographic experts group (JPEG) 
image compression standard algorithm, a motion picture 
experts group (MPEG) image compression standard 
algorithm, and a H.261 video compression standard algo- 
rithm. 

41. A discrete cosine transform-based digital image com- 
pression method having the step of generating a rate- 
distortion-optimal quantization table as claimed in claim 32 
and further comprising the steps of: 

(a) dividing input digital image data into fixed length 
image blocks and generating a corresponding block of 
discrete cosine transform coefficients for each fixed 
length image block by applying a discrete cosine trans- 
form to each fixed length image block: 
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(b) generating a corresponding block of quantized discrete 
cosine transform coefficients for each block of discrete 
cosine transform coefficients by quantizing each block 
of discrete cosine transform coefficients with the rate- 
distortion-optimal quantization table; and 

(c) entropy coding each block of quantized discrete cosine 
transform coefficients to form a corresponding com- 
pressed block such that a sequence of all of the com- 
pressed blocks forms the compressed digital image 
data. 

42. The method of claim 41 wherein the entropy coding 
step comprises entropy coding each block of quantized 
discrete cosine transform coefficients with Huffman codes to 
exploit similarities across the quantized discrete cosine 
transform coefficients. 

43. A method for operating a discrete cosine transform- 
based digital image compression and decompression system 
having the steps of discrete cosine transform-based digital 
image compressing as claimed in claim 41 and further 
comprising these steps of: 

(a) transmitting the compressed digital image data and a 
rate-distortion-optimal quantization table over a com- 
munication channel; 

(b) receiving the compressed digital image data and the 
rate-distortion-optimal quantization table from the 
communication channel; 


(c) decoding each entropy-coded compressed block in the 
received compressed digital image data into a corre- 
sponding block of quantized discrete cosine transform 

coefficients; 

5 

(d) generating a corresponding block of discrete cosine 
transform coefficients for each block of quantized dis- 
crete cosine transform coefficients based on the 
received rate-distortion-optimal quantization table; and 

10 (e) generating a corresponding decompressed block of 

fixed length digital image data for each block of 
discrete cosine transform coefficients by applying an 
inverse discrete cosine transform to each block of 
5 discrete cosine transform coefficients such that a 
sequence of all of the decompressed blocks forms 
decompressed digital image data. 

44. The method of claim 43 wherein the steps of trans- 
mitting and receiving utilize a communication channel 
20 selected from the group consisting of a serial communica- 
tion line, a parallel communication line, a data network, a 
public switched telephone network, an optically-encoded 
disk, a magnetically-encoded disk, and a radio communica- 
tion link. 

25 

* * * * * 



